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Abstract. We study non-linear structure formation in the presence of dark energy. 
The influence of dark energy on the growth of large-scale cosmological structures is 
exerted both through its background effect on the expansion rate, and through its 
perturbations as well. In order to compute the rate of formation of massive objects we 
employ the Spherical Collapse formalism, which we generalize to include fluids with 
pressure. We show that the resulting non-linear evolution equations are identical to 
the ones obtained in the Pseudo-Newtonian approach to cosmological perturbations, in 
the regime where an equation of state serves to describe both the background pressure 
relative to density, and the pressure perturbations relative to the density perturbations 
as well. We then consider a wide range of constant and time-dependent equations of 
state (including phantom models) parametrized in a standard way, and study their 
impact on the non-linear growth of structure. The main effect is the formation of 
dark energy structure associated with the dark matter halo: non-phantom equations 
of state induce the formation of a dark energy halo, damping the growth of structures; 
phantom models, on the other hand, generate dark energy voids, enhancing structure 
growth. Finally, we employ the Press-Schechter formalism to compute how dark energy 
affects the number of massive objects as a function of redshift (number counts.) 
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1. Introduction 

Observations of high-redshift SNIa imply that the expansion of the universe has been 
accelerating in the past few billions of years [H El [3l [4]. This is corroborated by 
at least three broadly independent observations: the angular spectrum of the cosmic 
microwave background temperature fluctuations O E], the galaxy-galaxy correlation 
function, which traces the spatial distribution of large-scale structure [3, [8], and the 
baryon acoustic oscillations [9]. Presently, the combined datasets favour a flat universe 
with Qrn — 0.27, where the remaining 73% of the energy budget is taken up by dark 
energy. 

These observations suggest that the dominant contribution to the present energy 
density of the Universe can be described by a dark {i.e., weakly or non-interacting) 
fluid with equation of state (EoS) Wde = Pde/Pde < —1/3. A particular case of such 
a substance would be the cosmological constant. A, for which w\ = —1. Many other 
models with Wdc 7^—1 have been proposed, usually in the framework of a scalar field 
("quintessence") or some other form of cosmic fiuid with negative pressure - see, e.g., 
|10| for a comprehensive review. 

A more direct approach to the phenomenology of dark energy has been recently 
adopted, in which the equation of state Wde is expressed in terms of a certain 
parametrization with respect to its time dependence [H], [121 1131 [H]- Although 
determining the equation of state as a function of redshift would probably not help 
to reveal the nature of dark energy, it could go a long way towards discriminating 
among existing models. Hence, one of the most important tasks ahead for observational 
cosmology is to gather sufficient data to successfully and inequivocally distinguish 
between this landscape of possibilities. As for theorists, the challenge is to determine 
in which additional ways dark energy may manifest itself in nature, apart from the 
acceleration of the overall expansion rate of the Universe. 

One of the ways in which dark energy changes the evolution of our local Universe 
is through its influence over the rates of formation and growth of collapsed structures 
(halos). Since all galaxies and quasars, as well as supernovae and putative sources 
of gamma-ray bursts, lie inside collapsed structures of some type or another, their 
distribution in size, space and in time will reflect to some extent the influence of dark 
energy. 

There are basically three mechanisms through which dark energy affects large-scale 
structure. First, the collapse of an overdense region due to gravitational instability 
is slowed down by the Hubble expansion drag, so any additional component which 
increases the expansion rate for the same value of the energy density will dampen the 
formation of collapsed structures. Second, as the accelerated expansion picks up speed, 
the large-scale gravitational potentials grow slower, then start to decay. This means 
that, as dark energy becomes the dominating dynamical component of the Universe, 
some large-scale overdensities will grow slower, and the process of gravitational collapse 
will even reverse itself at scales comparable to the Hubble horizon. And third, if dark 
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energy is not the cosmological constant then it must fluctuate both in time and in 
space. Hence, dark energy not only feels the gravitational pull of a matter halo, but 
it tends to form halos itself, thus influencing back those matter halos in a non-linear 
manner. Notice that the flrst two mechanisms affect collapsed structures only indirectly, 
through changes in the Hubble expansion rate, while the third mechanism depends on 
the clustering properties of dark energy. Since different models of dark energy can easily 
produce the same homogeneous expansion rate, but they hardly ever produce the same 
perturbations, the largest potential to probe the nature of dark energy possibly comes 
from such perturbative mechanisms. 

In a previous paper, two of the present authors have studied the influence at the 
background level (no dark energy fluctuations) of different parametrizations of the dark 
energy equation of state (EoS) in the evolution of dark matter perturbations and in the 
final number counts of dark matter halos p^. Our main purpose in this paper is to 
extend this analysis by studying non-linear structure formation including the possibility 
of dark energy fluctuations. 

Related approaches were recently developed by Nunes & Mota [16], Manera & Mota 
[T7] . Nunes, Silva & Aghanim |18] and Dutta & Maor [Hj, but those works considered 
scalar field dark energy. Here we focus instead on dark energy as described by some 
parametrization for the EoS as a function of redshift. This is more general than the 
scalar field approach, since the EoS is directly related to the physical observables most 
widely used to measure cosmic acceleration. Moreover, in contrast to [I6l[l7], we were 
able to investigate the non-linear regime of both dark matter and dark energy clustering 
consistently, and we have found that it has an important effect on the formation of 
massive objects (M > 10^^ Mq). Although the description of dark energy entirely in 
terms of its EoS may be too restrictive, the present work proves that the impact of 
dark energy perturbations on the formation of collapsed objects is both substantial and 
observable. 

This paper is organized as follows. In Sec. 2 we extend the traditional Spherical 
Collapse (SC) model, originally used only to describe gravitational collapse in the 
absence of pressure, to incorporate the possibility of coupled perturbations in 2 fiuids, 
namely pressureless dark matter and negative pressure dark energy. We verify this 
generalization showing under what conditions the SC model is equivalent to a Pseudo- 
Newtonian (PN) perturbation theory approach. In Sec. 3 we study the linear evolution 
of the generalized SC equations, and the effects of dark energy on the formation and on 
the initial stages of the evolution of matter halos. In Sec. 4 we analyze the fully non- 
linear system of SC equations, and how the formation of strongly non-linear (collapsed) 
matter halos both affects and is affected by dark energy halos. In Sec. 5 we show 
how our results can be included in a Press-Schechter formalism in order to derive the 
consequences to dark matter halos number counts. We conclude in Sec. 6. 
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2. The spherical collapse model and its generalizations 

The simplest (semi-) analytical tool to study non-linear structure formation is the SC 
model [201 . It has been shown that the SC equations can be actually derived from 
General Relativity, as long as shear does not play a significant role [2T] . 

Most studies about the impact of dark energy on structure formation were 
performed under the assumption that dark energy is uniformly distributed. In this 
case, where dark energy affects only background quantities, the SC model can be easily 
modified to incorporate dark energy effects. For instance, the abundance of rich clusters 
of galaxies estimated within the SC model was used to constrain the cosmological 
model and the properties of dark energy fluid in the context of the simplest case of 
a cosmological constant [22l [23], in the case of a constant w^e 7^—1 [SH [25], [26], [27l [28]. 
as well as the case of dynamical dark energy models with some parametrizations of 
waeit) m- 

However, the standard SC framework was originally designed to describe 
perturbations in pressureless matter, while we are interested in the effects of 
perturbations in an extra component whose pressure is very large and negative. If we 
want to study a gravitationally coupled system of matter (we do not distinguish between 
dark matter and baryons here) and dark energy, the SC model must be expanded beyond 
the realm of the Einstein-de Sitter model. 

Consequences of dark energy fluctuations in the studies of structure formation are 
more naturally incorporated by introducing a scalar field with a suitable potential to 
model the dark energy component, such as the quintessence field. In this approach, the 
authors of Refs. [16], HTl [29], [30] proposed an extension of the SC equations that take into 
account fiuctuations in the dark energy field for minimally and non-minimally coupled 
quintessence field. 

It is often more convenient, and completely equivalent at the background level, to 
introduce a time-dependent parametrization for the dark energy EoS, Wde{z). Since it is 
possible to reconstruct the scalar field potential from a general parametrization of dark 
energy, or directly from the EoS Wde{z) [Ml [32], the two approaches are in fact closely 
related. Our goal in this section is to generalize the SC model with a fiuid description 
of dark energy in order to include the possibility of dark energy fluctuations. 

We will check this generalized SC model with the results from the Pseudo- 
Newtonian (PN) approach to cosmology [33] for perfect fluids with pressure. The 
advantages of the PN framework are that it is both simpler than full-blown non-linear 
General Relativity (GR), and more intuitive. Crucially, it is in good agreement with 
GR in the linear regime [34l [35]. As we will see, the PN approach is also particularly 
useful if we want to keep contact with the description of dark energy in terms of a 
parametrization for its EoS, and it can be easily generalized to a multi-fluid system. 

The only remaining question is whether the two approaches agree with each other. 
Next we verify under which conditions the SC model is equivalent to the PN approach. 
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2.1. Spherical Collapse 

The continuity equation for a single perfect fluid j with background density pj and 
pressure pj = WjPj is given by: 

Pj + 3Hpj (1 + Wj) = , (1) 

where H = a/a is the Hubble parameter. Consider now a spherically symmetric region 
of radius r and with a homogeneous density pc^ (a top-hat distribution). Suppose that, 
at time t, Pcj{t) = Pjif) + 5pj- If ^Pj > this spherical region will eventually collapse 
from its own gravitational pull, otherwise it will expand faster than the average Hubble 
flow, generating what is known as a void. The evolution of such simplifled spherical 
regions can be described in close analogy with the continuity Eq. ([I]), but now with 

Pcj Pcj ■ 

Pc, + (1 + = , (2) 

where h = r/r denotes the local expansion rate inside the spherical region. Note 
that, in principle, we could have different equations of state inside and outside the 
spherical region, Wc^ 7^ wj. In fact, the difference between the local and the background 
equations of state 6wj = — Wj can be related to the fluid's effective speed of sound, 
4s j = hjl^Pj, through: 

^^i = T^^(Ceff,-^.)- (3) 
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Usually, Cgg is regarded as a free parameter - although, rigorously, in perturbation theory 
the only other free parameter is the true sound speed of inhomogeneities, [36j. Here 
the sound speed c^g is deflned as the ratio between two independent perturbative degrees 
of freedom, so not only it is gauge dependent, but it may depend also on the initial 
conditions for those perturbations. Therefore, Cgg stands as a proxy for the pressure 
perturbations. 

For simplicity, and in order to make contact with the PN equations, we will consider 
the case where the EoS is the same inside the collapsing sphere and in the background, 
so we take 6wj = and thus c^^j = Wj. This situation can be readily obtained in cases 
such slow-rolling scalar field. 

It should be noted that, in principle, there are instabilities in the growth of 
inhomogeneous perturbations whenever the sound speed becomes negative. However, 
within the spherical collapse model with a top-hat profile and the assumption of a 
space-independent Cgg, there are no pressure or density gradients, so no such problem 
of instabilities arises. 

By the same token as the first Friedmann equation, consider now the second 
Friedmann equation applied to the spherical region: 
r 4:7iG , , 

-^{Pc + ^Pc) ■ (4) 



r 
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Notice that the density and pressure that appear in Eq. ([4]) are the sum of densities and 
pressures of all contributing fluids, while the continuity Eqs. ([Il)-® are valid for each 
individual fluid (in the absence, of course, of direct couplings between those fluids.) 
It is useful to define the density contrast of a single fiuid species j by the relation: 

Sj + l = ^. (5) 

Differentiating this with respect to time we obtain: 

6j = 3{1 + 6j){H -h){l + Wj) , (6) 

where we assumed Wc^ = wj. Differentiating again with respect to time and employing 
the equations for the background and for the spherical region, we can derive the following 
non-linear evolution equation for 6j: 

Sj + (2H - 5j - A-kG (1 + Wj) (1 + 5j) J2 PkSk (1 + Swk) = 



4 + 3wj 
3(1 + Wi 
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1 + S. 



(7) 



Notice that we admit the possibility of a time-dependent EoS. For a system of n 
fluids, we must consider n equations such as ((71), all coupled gravitationally through 
the term proportional to Newton's constant. Although they are not derived rigorously 
from General Relativity, we will see next that these equations find support in the PN 
approximation to gravitational interactions. 

2.2. Pseudo-Newtonian Cosmology 
PN cosmological 

+ Vr ■ (Ujpj) + PjVr ■ Uj = , (8) 



Consider now the PN cosmological model, described by the equations [33] : 

dt 



^"^^ ' ^U^.Wr)uj = -Wr^-^^ , (9) 



dt \ ^ } ■' pj+ pj 

V^$ = 47rG'5^(pfc + 3pfc) , (10) 

k 

where pj, pj, Uj and $ denote, respectively, the density, pressure, velocity and the 
Newtonian gravitational potential of the cosmic fiuid. These equations are, respectively, 
generalizations of the continuity equation, of Euler's equation (both valid for each fiuid 
species j), and of Poisson's equation (which is valid for the sum of all fiuids.) 

Cosmological perturbations are introduced by admitting inhomogeneous deviations 
away from the background quantities: 

Pi = POj (i) + 5pj(^,i) , (11) 

Pj = POj{t) + Spj{x,t) , (12) 

Uj = MOj (t) + Vj{x, t) , (13) 

$ =$o(t)+0(^,t) • (14) 
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Changing to comoving coordinates, x = f/a, (henceforth V refers to gradient 
with respect to comoving coordinates x) and using 6j = Spj/po^, we find the following 
equations for the perturbed quantities: 



5j + 3H {clfij - Wj) Sj 



Vj + Hvj + 



[l + wj+{l + cls,)6j] 



V ■ Vi Vi ■ V6i 



V0 



-cSj 



,V5 



k 



(15) 
(16) 
(17) 



where cl^j = 6pj/6pj is the effective sound speed of each fiuid. In order to obtain these 
equations we have assumed that wj and Cgg^ are functions of time only. 

Notice that Eqs. (fT5ll -(fT7ll are valid even if 6j is not small, so we can use them to 
follow the evolution of a collapsing region well into the non-linear regime. In fact, the PN 
equations of motion become a better approximation as the size of the system shrinks due 
to gravitational collapse. This is easy to see by noticing that in most collapsed regions 
of the Universe the density contrast 6j may be extremely large, but the gravitational 
potentials are small, <^ 1, and the local relative (peculiar) velocities are almost never 
relativistic. Hence, the PN equations may be a poor approximation at the moment of 
turnaround (when a spherical region breaks away from the Hubble fiow) for the scales 
comparable to the Hubble horizon at the time of turnaround, but for all other scales 
and epochs it is a good approximation that becomes progressively better as the system 
collapses. 

In order to simplify the PN equations, it is useful to define: 



and 



/. = V ■ 



V0 



+ 



a (1 + Wj + 6j + clsj 
so that by taking the divergence of (fT6l) we obtain: 

9j + H9j + Cj = -fj . 
We also define: 



6,) 



1 



Wj + (l + C^sj 



(18) 
(19) 

(20) 

(21) 

(22) 
(23) 



and by neglecting the term vj ■ V6j, which is of order of w|/c^, we can cast Eq. (fTSj) in 
the form: 



L + A. + ^B. = . 



(24) 
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Taking the partial derivative of (l24l) with respect to time, using Eq. ( l2Ti l to eliminate 
9j and Eq. (l24l) to eliminate 9j we get: 

+ A, + (a, + 5,) (^2H - j - ^ (/. + C,) = 0. (25) 

In the appendix we explicitly show that, for a single fluid with clj-j- = w = const, PN 
and GR at the linear level differ only by a decaying mode. 

Now we try to make contact with the SC equations. In order to reproduce Eq. ((Tj) 
we must, first of all, assume that the velocity profile is consistent with the hypothesis of 
spherically symmetric collapse of a top-hat inhomogeneity, i.e., Vj = 9(t)/3x. Second, 
we also have to assume that Cg^^ = Wj. Notice that because the intrinsic non-adiabatic 
pressure Tj ~ 6pj — c^jSpj, where the adiabatic sound speed is c^ j = j-, our choice 
implies some amount of intrinsic entropy perturbations for the dark energy fiuid. 

With these choices we have: Aj = 0, Bj = {1 + wj) {1 + 6j) and Cj = 9"^ /3a. 
Notice that, with this velocity field, the LHS of equation (i2Ti ) is identical to that of 
Raychaudhuri's equation when we assume that 9j is a function of time only: 

+ + ^ = -h ■ (26) 

This equation reduces to the one found in [2T] if we neglect the gradients of the density 
contrast in fj. Since we are considering the spherical collapse of a top-hat distribution 
(which is homogeneous inside the radius r) , the terms including V5j which appear in fj 
vanish. Under these conditions, Eq. (I25ll reduces to the equation for SC, Eq. ([7]). 

It is interesting that in fact we were forced to assume both that = Wj in the SC 
formalism, and that Cg^^- = Wj in the PN formalism, in order that the two frameworks 
would result in identical equations. This is a further motivation for our choices of 
6wj = and c^gj = wf only in this scenario we can trust that the physics of non-linear 
spherically symmetric collapse is well described by our dynamical equations. In order 
to describe a more general situation probably neither approach is suited, and one would 
be forced to resort to full-blown General Relativity. However it is possible that the 
numerical differences between SC and PN for other choices of 5wj and c^^j are small. 

2.3. Equations for non-linear spherical collapse in the presence of dark energy 

We obtained non-linear differential equations that characterize the growth of spherically 
symmetric perturbations in fiuids with arbitrary time-dependent equations of state. 
These equations are coupled through the gravitational interactions. We saw that 
both the PN and SC approaches agree with General Relativity for a pressureless 
fiuid; furthermore, we have shown that they agree with each other in the case where 
^eff j — '^cj = Wj and if the density profile is a top-hat (V5j = .) 

Particularizing to a model with only non-relativistic matter and dark energy, in 
which the latter is characterized solely by its EoS, the top-hat spherical regions evolve 
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according to a system of equations equivalent to ([7]): 

L + 2HK - , = ^ (1 + 5m) [n^K + n^Jdc (1 + 3w7dc)] ,(27) 

3(1 + dm) 2 



1 + Wdc 



4 + 3Wde 

3(1 + Wdc) 



1+5, 



dc 



^ (1 + Wde) (1 + 5dc) i^mSm + ^deSde (1 + 3Wdc)] , (28) 

where 6^ is the density contrast in matter and 6de is the density contrast in the dark 
energy component. These are the equations we will study in the following sections. 

3. Solutions in linear regime 

The linear regime of cosmological perturbations is valid for all scales during the radiation 
era, and for most scales during the matter era up until very recently. The initial stages 
of the process of gravitational collapse are indeed very well described by the linear 
regime for all but the smallest scales. Since this is a simple system which can be studied 
almost entirely with analytical tools, it is useful to try and extract some physics from 
Eqs. (l27l) - (l28l) while they are still in the linear regime. Neglecting the 0{6'^) terms we 
obtain: 

K + 1H5m = [^^mf^m + ^deSde (1 + 3Wde)] , (29) 

5de+(2iJ--^^)5de (30) 



1 + Wd, 



= — ^ (1 + Wde) [^mSm + ^de^de (1 + 3Wdc)] ■ 

In principle, we can employ any given parametrization for dark energy as a function 
of time or redshift, but in order to find closed analytical formulas we initially take 
Wde =constant. We start by solving Eqs. (I29l)-(l30ll well inside the matter-dominated 
period {z = 10^), when it is a good approximation to assume that f2de ~ and f2m ~ 1. 
Changing the time variable to the scale factor a, the equations become: 

«o + |^-^(l + ™d,)i„ = 0, (32) 

where a prime denotes derivative with respect to a. 

As is widely known, in this case the solution for the matter density contrast is 
Sm{<^) = Ci a + 6*2 a"^''^, where Ci and C2 are arbitrary constants. Neglecting the 
decaying mode of the matter density contrast, Eq. ([32l) then has the solution: 

Sde = Ci (1 + Wde) a + C3 = (1 + Wde) K + C3 . (33) 

It is interesting to note that the adiabatic condition is 6dc = (1 + Wde) 5m- Hence, 
any value C3 7^ implies a non-adiabatic initial condition - i.e., in such a case 
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the perturbations have an isocurvature component. However, local non-adiabatic 
perturbations are unstable, since they correspond to pressure gradients between the 
internal and external parts of the spherical region. These pressure gradients must 
eventually cancel out, so any non-adiabatic component must decay during the evolution 
of the perturbations, leaving only the usual adiabatic (curvature) fluctuations. 

Notice that the condition of adiabaticity is different in the usual dark energy 
models (wde > —1) compared to phantom models {wde < —1): for phantom models, 
adiabatic initial conditions mean that any initial overdensity in matter is matched by 
an underdensity in dark energy, and vice-versa. So, for example, take a phantom dark 
energy model and a positive density perturbation in dark matter. If initially the dark 
energy perturbation is also positive, then the pressure gradients will cause the dark 
energy halo to decay, then to turn it into a void, thus switching the sign of the dark 
energy perturbation - see [29j for a similar switching effect. 

The effect of dark energy perturbations on the evolution of dark matter 
perturbations is easily understood from Eq. (I29l) : dark energy perturbations become 
a source for dark matter perturbation. Since (1 + 3wdc) < 0, a dark energy overdensity 
decreases dark matter clustering, which is intutitive since a local concentration of dark 
energy would speed up the acceleration in that region. The opposite holds for a region 
with a dark energy underdensity. We show some examples in the next subsections. For 
the remainder of the paper we adopt Qm = 0.25, Qde = 0.75 and ho = 0.72. 

3.1. Constant Wde : non-phantom models 

In the case of non-phantom dark energy we have 1 + Wdc > 0, and therefore Eq. (l33l ) 
implies that a region containing a matter overdensity {6^ > 0) induces a dark energy 
overdensity (^de > 0) in that same region. Conversely, a matter underdensity region 
(5m < 0) induces a dark energy underdensity {6de < 0). Hence, in the non-phantom case 
a halo of dark matter induces a halo of dark energy, and a void of dark matter induces a 
void of dark energy. This behaviour is generic if we limit the scope of initial conditions 
to adiabatic perturbations - as predicted by inflation and conflrmed by WMAP [6]. 

In order to study the impact of dark energy fluctuations on the growth of dark 
matter perturbations we show in Fig.[I]the evolution of 5m with and without the inclusion 
of dark energy fluctuations, for two different values of the EoS: Wde = —0.9 and —0.8. 
We use adiabatic initial conditions at Zi = 1000, with S'^i^Zi) = —6m{zi)/{l + Zi). 

The initial condition on the derivative of the density perturbation comes from the 
assumption that dark energy is negligible initially. As the figure shows, the inclusion of 
dark energy perturbations actually suppresses the growth of dark matter perturbations 
in this case. The differences in the linear regime are roughly proportional to 1 + Wde- 
In this same figure we also show the growth of the dark energy perturbation, which is 
much smaller than the matter fiuctuation, as expected, but also tends to form a dark 
energy halo. 
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Figure 1. Growth of matter perturbation with (solid lines) and without (dashed lines) 
dark energy perturbations together with growth of dark energy perturbations (dotted lines) 
for Wdc = —0.8 (left panel) and Wdc — 0.9 (right panel). 



3.2. Constant w^e'- phantom models 

In the case of phantom dark energy we have 1 + w^e < 0, and therefore Eq. ( l33l l implies 
that a matter overdensity region (^m > 0), which will later become a dark matter halo, 
induces a dark energy density void (5de < 0), and vice- versa. Again, this behaviour is 
generic for purely adiabatic initial conditions. 

In Fig. [2] we show the effects of dark energy fluctuations on the growth of dark 
matter perturbations for two different values of the EoS, w^e = —1-1 and —1.2. As in 
the case of non-phantom dark energy, the differences are small and increase with larger 
values of |l + Wde|- However, contrary to the non-phantom case, fluctuations in phantom 
dark energy enhance the growth of dark matter perturbations. We also show the growth 
of the dark energy perturbation which, as expected, tends to form a dark energy void. 

3. 3. Varying Wde 

In the framework of single scalar held descriptions of dark energy it is impossible for 
the EoS to cross the so-called phantom barrier at Wde = —1 |37]. However, in our 
phenomenological approach we could in principle have a time- varying parametrization of 
Wde{z) crossing the phantom barrier. In fact, this is the case with many parametrizations 
adjusted to flt SNIa data [T4] . 

The existence of a phantom barrier is hinted in our approach by the presence of the 
term Wde (1 + u)de)~'^ in Eq. f l30l ). Although the divergence at Wde = — 1 is not necessarily 
fatal for the solutions of the differential equations, here we consider only dark energy 
parametrizations that are phantom or non-phantom during all times. 

We will study a parametrization of the dark energy EoS of the form [38] : 

Wde = Wo + wi{l - a) = Wo + wi— — , (34) 
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Figure 2. Growth of matter perturbation with (solid lines) and without (dashed lines) 
dark energy perturbations together with growth of dark energy perturbations (dotted lines) 
for Wdc = —1-1 (left panel) and Wdc — —1.2 (right panel). 



and we choose parameters Wq and Wi which are consistent with the 2a regions which 
are jointly constrained by observations of the CMB, supernovas and baryon oscillations 
- see, for instance, [39l [40] . 

In Fig. [3] we show the impact of dark energy flucutations for a variable EoS for both 
phantom and non-phantom cases. The results are similar to the constant EoS cases. 




Figure 3. Evolution of {z) in linear regime with dark energy adiabatic IC with 
Wo — —0.75 and wi = 0.4 (non-phantom, left figure) and with wq = —1.1 and wi = —1 
(phantom, right figure) including dark energy perturbations (full line) and without dark 
energy perturbations (dashed line). The growth of dark energy perturbations (dotted line) 
is also included. 
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4. Non-linear regime 

In the non-linear regime, as in the linear regime, we consider again only models 
that are phantom or non-phantom at all times. We solve Eqs. (I27ll -(l28l). and for 
brevity's sake we limit our scope to the parametrizations discussed in subsection I3.3[ 
{wo,wi) = (—0.75,0.4) and (—1.1,-1). In particular, we are interested in the impact 
of the dark energy fluctuations on the collapse of dark matter structures. We will see 
that the effects found in the linear case (of the order of a few percent) are amplified by 
the non-linear evolution. 

In Fig. [4] (left panel) we show the dark matter density contrast for initial conditions 
chosen such that the collapse (indicated by the divergence of the density contrast) of a 
spherical dark matter structure happens at the present time {z = 0). Physically, this 
would correspond to the formation of an object such as a supercluster. Dark energy 
fluctuations have a dramatic effect in this case: in the non-phantom case, that structure 
would have collapsed much earlier if the dark energy fluctuations had not been taken 
into account. In the right panel of Fig. [4] we also show a phantom parametrization, 
where initial conditions are chosen such that the collapse of the dark matter structure 
takes place today without dark energy perturbations. In this case, the inclusion of 
dark energy perturbations enhances the clustering of dark matter and cause that same 
structure to collapse earlier. 




Figure 4. Non-linear evolution of Sm (z) for dark matter with (solid) and without (dashed) 
dark energy perturbation for the case of a non-phantom (left panel) and a phantom (right 
panel) parametrization. In the left panel, initial conditions are chosen such that dark matter 
with dark energy perturbations collapses today (z — 0), whereas in the right panel initial 
conditions are chosen such that dark matter without dark energy perturbations collapses 
today. 

In Fig. [5] we contrast the evolution of the non-linear densities in dark matter and 
dark energy against the evolution of the linearized densities. In that flgure the initial 
conditions are chosen such that the non-linear dark matter perturbation diverges at 
z = 0. Hence, the value of the linearized dark matter perturbation at = in this case 
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corresponds to the definition of the critical density contrast for that redshift, Sc{z = 0). 
In the next Section we will employ the critical density contrast as a function of redshift 
in order to to compute the Press-Schechter function. 




Figure 5. Non-linear evolution of 5 (z) for dark matter (solid line) and dark energy 
(dot-dashed line). The initial conditions are chosen such that dark matter collapses 
today (z — 0). Also shown are the linearized evolution of dark matter (dashed line) 
and dark energy (dotted line) perturbations. The left and right panels correspond to the 
parametrizations {wo,wi) = (—0.75,0.4) (left panel) and (—1.1,-1) (right panel.) 



5. Number Counts of Dark Matter Halos 



We will use the Press-Schechter approach [4T] in order to estimate the number 
counts of dark matter halos for different bins of redshifts and halo masses. More realistic 
mass functions could be used (see, for instance, [42]) but our intention here is simply 
to point out how number counts differ in scenarios with dark energy compared to the 
standard ACDM model. We believe these differences would be essentially the same 
in more sophisticated models of nonlinear structure formation, such as the model of 
ellipsoidal collapse. 

The Press-Schechter formalism assumes that the fraction of mass in the universe 
contained in gravitationally bound systems with masses greater than M is given by the 
fraction of space where the linearly evolved density contrast exceeds a threshold 6c, and 
that the density contrast is normally distributed with zero mean and variance a^{M) 
- the root-mean-square value of the density contrast 6 at scales containing a mass M. 
Therefore, it is assumed that for a massive sphere to undergo gravitational collapse at 
a redshift z, its linear overdensity should exceed a certain threshold 6c, defined as the 
linearly evolved density contrast at the instant when the non-linear density contrast 
associated with the mass M diverges(«. e., at the moment of collapse.) Since each scale 
M collapses at a given redshift (bigger masses, corresponding to larger scales, collapse 
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later), the critical density contrast is a function of redshift, 6c = Sc{z). Nevertheless, 
notice that only linear quantities are used in this formalism. For a review of the Press- 
Schechter formalism see, for instance, [43]. 

In order to illustrate how dark energy affects gravitational collapse, in Fig. [6] we 
show 6c{z) in a few dark energy scenarios. Notice that the inclusion of dark energy 
perturbations has a dramatic effect in Sc{z), with a substantial suppression in the non- 
phantom case and a substantial enhancement in the phantom case. If dark energy 
perturbations are not included we reproduce the results of [T5] . 
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Figure 6. Values of the linear critical density contrast 8c as a function of the collapse 
redshift. The solid line is the usual ACDM result. The values which do not include dark 
energy perturbations are shown for phantom (dashed line) and non-phantom (dot-dashed 
line) cases. The values with dark energy perturbations are shown for phantom (long 
dashed line) and non-phantom (long dot-dashed line). The non-phantom and phantom 
cases correspond to the same parametrizations used in Fig.[3l 



These assumptions lead to the well-known analytical formula for the comoving 
number density of collapsed halos of mass in the range M and M -|- dM at a given 
redshift z\ 



:{M,z) = -J exp 



2a2(M,z) 



(35) 



dM^ ' ' V TT M a{M,z) dM 

where is the present matter density of the universe and 6c{z) is the linearly 
extrapolated density threshold above which structures collapse, i.e., 5c{z) = 6ii^{z = 

Zcol)- 

The quantity: 

a{M,z) = D{z)aM (36) 

is the linear theory rms density fluctuation in spheres of comoving radius R containing 
the mass M, where D{z) = 6m{z)/Sm{z = 0) is the linear growth function obtained 
from Eq. (1291) . The smoothing scale R is often specified by the mass within the volume 
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defined by the window function at the present time, see e.g. [44]. In our analysis we 
use the fit given by 



j^^X-7(M)/3 

^M = CTs(—) , (37) 



MsJ 

where Mg = 6x lO^^^l^^^ Mq is the mass inside a sphere of radius Rg = 8/;.~^Mpc, and 
as is the variance of the over-density field smoothed on a scale of size Rs. The index 7 
is a function of the mass scale and the so-called shape parameter, T = Q!f)h Q-^b-^b/^ 
{VLh = 0.05 is the baryonic density parameter) [23j 



(0) 
AI 



7(M) = (0.3 r + 0.2) 



2.92 + 1 log 

3 ^\MsJ 



(38) 



For a fixed erg (power spectrum normalization) the predicted number density of 
dark matter halos given by the above formula is uniquely affected by the dark energy 
models through the ratio 6c{z)/D{z). In order to compare the different models, we will 
normalize to mass function to the same value today, that is, we will require: 

_ ScModjz = 0) 

(^8,Mod — —J ■p^<^8,A^ I'jyj 

Oc,a{^ = 0) 

where the label Mod indicates a given model and we use ag^A = 0.76 [6]. 

The effect of dark energy on the number of dark matter halos is studied by 
computing two quantities. The first one is the all sky number of halos per unit of 
redshift, in a given mass bin: 

.rhm dN f ,^ /"^^^-p dn dV 



dz j^^ J Mini d^ dzdfl 



where the comoving volume element is given by: 

r\z)/H{z), (41) 



dV{z) 2, 



dzdfl 

where r{z) = H~^{x)dx is the comoving distance. Note that the comoving volume 
element depends only on the cosmological background and is identical for models with 
and without perturbations in dark energy. The diferences in the number counts when 
one includes dark energy clumping are due to modifications in 6c{z) and D{z). 

The second quantity we compute is the all sky integrated number counts above a 
given mass threshold, Minf, and up to redshift z [T6] : 

^(^■^'>^^^//n£f^^<^A«.'. (42) 

Our knowledge of both these quantities for galaxy clusters will improve enormously 
with upcoming cluster surveys operating at different wavebands, such as the South Pole 
Telescope [45] . 

We can now examine the modifications caused by a clustering dark energy 
component on the number of dark matter halos with the same observable computed 
in the standard ACDM model. First we show how the different equations of state 
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3x10 




Figure 7. Evolution of number counts in mass bins with redshift for objects with masses 
within the range lO^^ < M/{h-^MQ) < lO" (top left panel), lO" < M/{h~^MQ) < 
IQi^ (top right panel) and 10^^ < M/{h-^MQ) < 10^^ (bottom panel). The solid 
line corresponds to the fiducial ACDM result. The number counts including dark energy 
perturbations are shown for non-phantom (long dot-dashed line) and phantom (long dashed 
line) models. The results without the inclusion of dark energy perturbations are also shown 
for non-phantom (dot-dashed line) and phantom (dashed line) models. 

impact the number of dark matter halos in given mass bins [Minf , Mgup] typical of the 
present-day cosmological structures, namely [10^^, 1O^'']/;,~^M0, [10^'', lO^^]/i~^M0 and 
[10^^,10^^] h~^MQ. The number counts in mass bins, A/"*^™ = dN/dz, obtained from 
(l40l) . are shown in Fig. [7l In each panel we plot the actual number counts together with 
the number counts computed for a fiducial ACDM model (solid lines), for each mass 
bin. Notice that the more massive structures are less abundant and form at later times, 
as it should be in the hierarchical model of structure formation. There is a slight shift of 
the peak redshift for structure formation in the distinct dark energy models considered. 
The differences with respect to the ACDM model become more significant in the bins 
with larger masses - but, of course, given the small number of such massive objects, the 
uncertainty due to shot noise also becomes increasingly important. 

Another important observable quantity is the integrated number of collapsed 
structures above a given mass, Eq. (I42l) . We present results for the integrated number 
counts of structures with masses above 10^^h~^MQ, IO^'^H^^Mq, and IO^^H^^Mq (we 
cut-off the integration at Mgup = lO^^/i^^M©, as such gigantic structures could not in 
practice be resolved today.) The results are displayed in Fig. [H always compared with 
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Figure 8. Evolution of the integrated number counts up to redshift z for objects 
with M > W^^h^^MQ (top left panel), M > lO^/i-^M© (top right panel) and 
M > 10^^ft~^Mo (bottom panel). The lines correspond to the same cases as in Fig. [71 



the results for the fiducial ACDM model (solid lines.) Notice that the integrated number 
has a plateau that refiects the epoch of structure formation for a given mass. In other 
words, there is no formation of structures with mass above 10^^ H'^Mq, IO^^H'^Mq, and 
10^^h~^MQ for redshifts roughly above z = 2, 1.5 and 0.7, respectively. Again we find 
large differences compared to the ACDM model when dark energy perturbations are 
included. 

6. Conclusions 

Our main goal in this paper was to study the effects of including dark energy 
perturbations in the evolution of matter perturbations in the linear and in the non-linear 
regimes. Since we do not know what dark energy really is, we developed a formalism 
whereby we can directly use a parametrization of its equation of state in order to address 
this issue. 

We have shown that the spherical collapse and the pseudo-newtonian approaches 
to the study the non-linear evolution of dark matter and dark energy perturbations are 
equivalent when one adopts an effective speed of sound Cgg = w. In the languange of 
spherical collapse, this is equivalent to assuming that the equation of state is the same 
inside and outside the collapsed region. 

We found distinct behaviours in the evolution of the dark matter perturbations 
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for phantom and non-phantom forms of dark energy. Inclusion of dark energy 
perturbations inhibits the growth of dark matter perturbations for the non-phantom 
case but it enhances this growth in the phantom case. The reason is that dark matter 
overdensities lead to dark energy overdensities in the non-phantom case, but they lead 
to underdensities in the phantom case. Due to its gravitationally repulsive nature, dark 
energy overdensities inhibit, while dark energy underdensities help, the growth of dark 
matter perturbations. 

This effect is small in the linear regime but becomes dramatic when studying the 
collapsed regions that have formed more recently. In particular, we found a large 
modification in the critical density 6c{z) even for moderatly low redshifts. 

We used the Press-Schechter formalism to estimate the modifications due to dark 
energy perturbations in observational quantities such as number counts of galaxy 
clusters, which refiect the formation and distribution of dark matter halos. We found 
that there are large deviations compared to the standard ACDM model, which are more 
significant for the larger structures. We expect that these large deviations are a general 
consequence of taking into account the non-linear dark energy perturbations that are 
gravitationally coupled to the dark matter perturbations. 

Although our use of the EoS to describe dark energy perturbations clearly 
constitutes a particular case, our choice was guided by the equivalence between the 
SC and PN approaches. However, we believe that more general models can still be 
described consistently in both approaches |46j. Hopefully future data on number counts 
of galaxy clusters will be able to discriminate among different models of dark energy. 
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Appendix A. 



In this paper we have analyzed the non-linear evolution of dark energy and dark matter 
perturbations using the PN. One can ask about the validity of such a model when 
pressure is taken into account. Here we explicitly show that, up to linear order, pseudo- 
newtonian equations are in good agreement with general relativity perturbation theory. 

The relativistic equation for the gauge invariant perturbation, for a single perfect 
fiuid can be written as l47l: 



5+2H 



w — — 



54 



{3w^ -8w-l + 6c1) S = — ^clsS .(A.l) 



Neglecting the term on RHS, which can be understood as a large scale approximation 
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or a top-hat profile, and taking w = const = cj. we have 

3w 



6 + 2H 



1 - 



ott2 

5 + — - {3w^ -2w-l)6 = 0, (A.2) 



nsingH = 2/3{l + w)t: 



V 4 — 6w 6 6w^ — Aw — 2 6 ^ , . , 

The solution to the relativistic equation above can be written as: 

5 = C'^t2(l+3«;)/3{l+t«) _j_ Q^^-l+2w/l+w _ 

In PN, Eq. (l25l) . up to linear terms, with c^g = ti? = const and 5pAr with no spatial 
dependence, is written as: 

ott2 

6pN + 2H6pN - — - (1 + w){l + 3w) SpN = 0. (A.5) 
and its solution is: 

The relativistic and pseudo-newtonian solutions differ only in the decaying mode. 
Note that for dust, w = 0, the two solutions coincide, even for the SC model [2X11481 [49] . 
Although both theories seem to be in good agreement, we must be watchful of phantom 
models. In this case the usual relativistic decaying mode ~ f--i+2w/i+w bg^omes a 
growing mode, while in PN it is always decaying. This behaviour is expected when 
phantom dark energy starts to dominate and shows that growing perturbations phantom 
dark energy increase the gravitational potencial, unlike one should expect from a 
homogeneous phantom dark energy model. For all other values of — 1 < w < 1 the 
relativistic solution ~ f-Ww/i+w -g ^ decaying mode. Hence, if the sub-dominant mode 
is not irrelevant for some reason, then the PN may not be reliable for phantom models at 
times when dark energy is strongly dominating the background evolution, i.e, fide ~ 1- 
Fortunately, this situation does not arise in our study. 
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